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无 单元 Galerkin 方法 中 周期 边界 条 件 的 处 理 * 


王晓东 ， 欧阳 洁 ， 苏 进 
(西北 工业 大 学 应 用 数学 系 ， 西 安 710129) 
摘 E 本 文 研究 了 无 单元 Galerkin 方法 中 周期 边界 条 件 的 处 理 技术 ， 将 Lagrange 乘 子 法 用 于 周期 边界 


条 件 的 处 理 . 数值 计算 结果 表明 ， 该 方法 具有 较 高 的 计算 精度 . 另外， 它 与 无 单元 Galerkin 方 
法 中 本 质 边界 条 件 处 理 的 Lagrange 乘 子 法 具有 统一 性 ， 对 于 周期 、 本 质 混合 型 边界 条 件 的 处 理 


尤为 方便 . 
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1 引言 


一 些微 分 方程 的 解 在 整个 求解 区 域内 会 具有 某 种 周期 性 ， 例 如 : 电磁 场 问 题 、 某 些 流动 问 
题 、 聚 合 物 分 子 的 取向 分 布 问题 等 .针对 这 类 问题 ， 我 们 可 以 利用 周期 性 将 问题 限定 在 一 个 周 
期 内 进行 求解 ， 从 而 降低 问题 规模 并 大 大 简化 计算 . 如 何在 限定 区 域 上 给 定 边界 条 件 ， 将 直接 
关系 到 微分 方程 的 求解 . 通常 可 以 给 定 两 种 边界 条 件 ， 一 种 是 自由 边界 条 件 ， 另 一 种 是 周期 边 
界 条 件 . 当 解 同时 具有 某 种 对 称 性 ， 并 且 可 以 确定 对 称 线 或 对 称 面 的 位 置 时 ， 可 以 将 求解 区 域 
取 为 两 个 相 邻 对 称 线 或 对 称 面 之 间 的 区 域 ， 并 且 由 于 解 在 对 称 线 或 对 称 面 上 导数 为 零 ， 所 以 可 
以 给 定 自然 边界 条 件 . 而 当 解 不 具有 对 称 性 或 无 法 确定 对 称 线 或 对 称 面 的 位 置 时 ， 则 一 般 只 能 
给 定 周期 边界 条 件 . 

无 单元 Galerkin (element free galerkin, EFG) 方法 [1- 引 作为 一 种 新 型 的 微分 方程 求解 方法 ， 
近年 来 受到 很 多 学 者 的 关注 .该 方法 不 依赖 于 网 格 或 单元 ， 具 有 易于 进行 自 适应 分 析 、 容 易 构 
造 高 阶 可 导 的 近似 函数 、 计 算 精 度 高 等 众多 优点 . 由 于 存在 这 些 优 点 ，EFG 方法 非常 适合 求 
解 诸 如 复杂 流体 多 尺度 模拟 此 6 中 具有 周期 性 边界 条 件 的 Fokker-Planck 方程 10] 等 问题 . 它 
的 高 精度 可 以 很 好 地 避免 Fokker-Planck 方程 求解 过 程 中 由 于 误差 积累 导致 的 取向 分 布 概率 密 
度 函 数 失 去 守恒 性 的 问题 . 要 实现 EFG 方法 在 这 类 问题 中 的 应 用 ， 首 先 就 必须 要 解决 周期 边 
界 条 件 的 处 理 问 题 . 

EFG 方 法 使 用 移动 最 小 二 乘法 3 构造 形 函 数 ， 而 移动 最 小 二 乘法 不 具有 Kronecker 
delta 函数 的 性 质 ， 所 以 有 限 元 等 方法 中 周期 边界 条 件 的 处 理 方 法 不 能 直接 用 于 EFG 方法 . 本 
文 借鉴 在 EFG 方法 中 使 用 Lagrange 乘 子 法 Qt13] 处 理 本 质 边 界 条 件 的 思想 ， 提 出 了 周期 边界 
条 件 处 理 的 Lagrange 乘 子 法 .该 方法 精度 高 ， 并 且 与 本 质 边 界 条 件 处 理 的 Lagrange RFR 
有 统一 性 ， 易 于 实现 . 
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* 基 金 项 目 : 国家 自然 科学 基金 (10871159; 10590353); 国家 重点 基础 研究 发 展 计 划 (973) 项 目 (2005CB321704). 
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2 EFG ARAN BRIE 


EFG 方法 使 用 移动 最 小 二 乘法 进行 未 知 函 数 的 逼近 和 检验 函数 的 构造 移动 最 小 二 乘 
法 03 (moving least squares, MLS) 是 局 部 近似 与 最 小 二 乘 相 结合 的 产物 ， 该 方法 直接 基于 空 
间 离 散 点 来 进行 近似 ， 不 需要 对 区 域 进行 网 格 划分 . 

在 区 域 Q 中 ， 未 知 标量 函数 w(z) 的 移动 最 小 二 乘 逼 近 可 以 表示 为 


u(x) = >》 Dj(x)aj(x) = p” (x)a(x), (1) 
j=l 
其 中 及 为 节点 间 的 平均 距离 ，pT(x) 为 基 函 数 向 量 ，m 为 基 函 数 个 数 . 二 维 空间 中 ，x = 
[x1, zz] 工 ， 基 函数 pT7(x) 可 以 取 为 


常数 基 pT = (1)， 
线性 基 p? = (1, £1, £2), 
二 次 基 p7 == (1, £1, £2, £7, 2122, 23), 


R (1) F a(x) = {a1(x), a2(x),--- ,am(x)}7 是 与 x 有 关 的 未 知 系数 向 量 ， 为 了 确定 a(x)， 定 义 
上 述 允 近 误差 的 加 权 Lo 范 数 为 


N 
R(x) = So w(x = x) [p7 (xi)a(x) — ui)’, (2) 
{=L 
其 中 为 节点 总 数 ，uwi = u(x) 为 函数 在 节点 xi RAE, wa- xi) 为 节点 权 函 数 ， 一 般 
取 为 三 次 样 条 函数 ， 即 


2 — Ar? + 4r3, ee 
w(x—x;)=w(r)= 4 4—4r44r?9- 40, Ler, 
0, r>1. 
为 保证 ut (x) 为 最 佳 近似 ， 需 使 R(x) 达到 最 小 ， 即 满足 
OF L A(x)a(x) - B(x)U = 0， (3) 
这 里 
N 
A(x) = 2, w(x 一 xi)p(xi)p (xi), 
B(x) = [w(x — x1)p(x1), w(x — x2)p(X2),.… ,w(x 一 xN)p(xN)]， 
U = {uu uy}. 
从 (3) 式 中 求 得 a(x) 并 代入 (1) 式 中 ， 便 可 以 得 到 函数 w(x) 的 MLS WEA 
u(x) = p?(x)A7*(x)B(x)U = ®7(x)U, (4) 
其 中 


T(x) = [px(z),paz(z)… ,pn(z)]. 
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3 EFG 方法 的 方程 离散 


假设 £ 为 定义 在 区 域 人 9 上 的 任 一 微分 算 子 ，% 为 满足 一 定 光滑 性 要 求 的 未 知 函数 ，f 为 已 知 
函数 ， 微 分 方程 的 一 般 形 式 可 以 写 为 





L(u) = f, (5) 
由 于 (5) 式 在 区 域内 任意 点 都 满足 ， 因 此 对 任意 函数 都 有 
(L(u), vo = (f, vo, (6) 


其 中 ， 函 数 v 称 为 检验 函数 ，(.,)89 AKRO 上 通过 函数 内 积 定义 的 泛 函 ， 即 
(ahja = | ghao, 
A (6) 称 为 微分 方程 (5) 的 等 效 积 分 形式 . 
一 般 地 ， 对 于 给 定 的 定 解 条 件 ， 方 程 (5) 很 难 精确 求解 ， 只 能 采用 数值 的 办 法 近似 求解 . 
BE ee! (x) 为 满足 定 解 条 件 的 方程 (5) 的 一 个 近似 解 ， 称 为 试探 函数 ， 它 可 以 表示 为 一 组 已 知 函 
数 (pi (x) 的 线性 组 合 ， 即 


N 
u(x) © w(x) = Y gis, (7) 
i=1 


ui = U(Xi) 为 待定 参数 .近似 解 u*(x) 一 般 不 能 精确 满足 微分 方程 (5)， 试 探 函 数 的 项 数 N 越 
多 ， 近 似 解 的 精度 就 越 高 ， 当 项 数 N 趋 于 无 穷 大 时 ， 近 似 解 将 收敛 于 精确 解 上 3l. 

如 果 对 任意 检验 函数 w，(6) 式 都 成 立 ， 此 时 近似 解 便 为 精确 解 ， 实际 上 ， 不 可 能 也 不 需要 
无 穷 多 个 检验 函数 ， 而 是 一 般 将 检验 函数 取 为 一 组 基 函 数 的 线性 组 合 ， 即 


M 
v= 》 cidi(x), (8) 
i=l 
其 中 M >N. 将 (7),(8) 两 式 代 入 (6) 式 ， 并 考虑 到 系数 ci 的 任意 性 ， 得 
N 
(£( > valu), o) = (f, 9;(x))o j = 1,2, „M, (9) 
i=1 
上 式 的 意义 是 通过 选择 合适 的 待定 参数 w:， 使 得 微分 方程 (5) 在 某 种 平均 意义 上 满足 . 
对 EFG IE, p(x) 通过 MLS 逼近 得 到 ， 试 探 函数 与 检验 函数 取 自 同一 函数 空间 ， 且 M = 


N, BẸ 
$3(x) = p;(x), j= 1,2,---,N, 


此 时 (9) 式 具体 为 
N 
(c( Daou) o) = (f, pi(x))o, 了 = 1,2,--- sN; (10) 
或 写成 本 
2, 〈《C(es(x))， pi(X)) ou = (f, £5(X))o> j=1,2,--- IN, (11) 


R (11) ÆA u 为 未 知 量 的 N 阶 代数 方程 组 . 

计算 时 可 以 使 用 分 部 积分 或 Green 公式 降低 £ 中 所 含 导数 的 阶 次 ， 从 而 降低 对 函数 wu*(x) 的 
连续 性 要 求 . 通过 定 解 条 件 对 代数 方程 组 (11) 做 相应 修正 后 求 得 w;， 将 wi 代入 (7) 式 即 可 得 到 
LAL u” (x). 
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4 ”周期 边界 条 件 的 处 理 
假设 微分 方程 (5) 具有 如 下 周期 边界 条 件 


u(T1) = u(L2), (12) 





其 中 边界 Ti 和 Ts 空间 位 置 不 同 ， 但 几何 形状 和 几何 尺寸 相同 ， 例 如 : 两 个 单 点 、 两 条 长 度 相 
同 的 线段 、 两 个 大 小 和 形状 均 相 同 的 平面 等 ，Ti 与 Fa 上 的 点 一 一 对 应 ，(12) 式 表 示 对 应 点 处 
未 知 函数 取 值 相同 . 由 于 (12) 式 在 边界 Ti 和 Ts 的 所 有 对 应 点 处 满足 ， 并 且 边 界 工 和 Ts 的 几 
何 形状 和 几何 尺寸 相同 ， 因 此 对 任意 的 入 都 有 如 下 的 积分 形式 成 立 


人 Auds 一 f Auds = 0, (13) 
Ta rs 
式 中 入 为 Lagrange 乘 子 [11,12]. 


FAD He BR BCA AS ST PRB u EES uw 时， 微分 方程 对 应 的 等 效 积分 形式 (也 称 Galerkin 变 分 
形式 ) 可 以 写 为 
61I(v) =f {duL(u) — fou}dQ = 0, (14) 
2 


将 (13) 式 以 变 分 的 形式 引入 到 (14) 式 中 得 到 


611, (u, A) = | tuc- foujao+ f ô\udsı 


+f Aĝudsı -f druds2 -f Aðudsz = 0, (15) 
Ti Te T2 


6u 和 6 入 分别 为 ww 和 入 对 应 的 变 分 函数 ， 且 bu e HHN), 5 € H°(Q), HP HQ), HO(Q) 为 1 
阶 和 0 阶 的 Sobolev 空间 . 对 任意 的 du Al OA, ME (15) 式 恒 成 立 的 入 必然 同时 满足 微分 方程 (5) 
和 周期 边界 条 件 (12). 另外 ， 由 于 (15) 式 等 价 于 


| {duL(u) 一 fou}an+ f Aðudsı -f Adudse = 0, (16) 
Q rı T2 

/ l 6 和 Xudsl 一 | duds = 0, (17) 
所 以 同时 满足 (16) 和 (17) RAY u BEBO AE. 


计算 时 ，w、 入 分 别 采用 MLS iB Al Lagrange 播 值 近 似 099， 并 且 选 取 bu 和 5 入 分别 为 MLS 
形 函 数 yi (x) M Lagrange 插值 基 函 数 N(x). 对 应 的 代数 方程 组 形式 为 


N 
> (leil), 25) qui 
i=l 
Np 


+ >， ( (Ne(x)， oOj(x))r ra) Mk = (f,93(*)) os j = 1,2,--- N, (18) 


k=1 


N 
` (Llpi(x)), Ne(X))p, r, u: =0, &=1,2,-+-, Np, (19) 
i=1 
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AP No 为 边界 Ti 或 Ta 上 的 节点 数 . (18) F (19) RÆ A u; AA, 为 未 知 量 的 N + N, 阶 代 数 方 
程 组 . ORG u; 便 可 得 到 满足 周期 边界 条 件 (12) 的 微分 方程 (5) 的 数值 解 . 

在 (15) 式 中 ， 周 期 边界 条 件 是 以 弱 形 式 引入 的 ， 即 周期 边界 条 件 是 在 平均 意义 下 满足 的 . 
也 可 以 令 边 界 Ti 51r 上 所 有 对 应 节点 满足 周期 边界 条 件 (12)， 即 


u”(xi) = u? (xi2), Xi ET1, Xi2 Er 1=1,2,.… , Np; (20) 


在 计算 过 程 中 ， 引 入 (20) 式 等 价 于 在 计算 (15) 式 中 的 边界 积分 时 采用 节点 积分 . 


5 ”数值 算 例 


5.1 一 维 情形 
对 均 质 棒状 分 子 聚 合 物 稀 溶液 ， 把 分 子 的 旋转 限定 在 平面 内 时 ， 聚 合 物 分子 在 微观 的 取向 
概率 分 布 函数 将 满足 如 下 的 一 维 Fokker-Planck 方程 [10] 


OY(6,t) _ 
Ot 





Wool, t) — 6g(C(n, 0)w(0,t)), (21) 


给 定 周期 边界 条 件 
Y0, t) = o(2n,t), (22) 


RP (0, HARADA PB, 0 € [0,27] 为 构 型 变量 ，« 为 宏观 流 场 速度 梯度 张 
量 ，De 为 无 量 岗 的 Deborah 数 ，C(i,0) =u x «u Hu = (cos(9),sin(9))7. 
在 二 维 前 切 流 场 和 单 轴 拉 伸 流 场 两 种 情形 下 ，«k 分别 取 


(o a Ps (3 0 ) 

0 0 0 =1 

计算 中 Deborah 数 取 De = 10, HER Y(6,0) = 1/2r， 取 向 空间 内 均匀 布置 101 个 节点 ， 每 个 
单元 上 采用 4 点 Gauss 积分 ， 时 间 步 长 取 为 0.005. Fokker-Planck 方 程 的 稳 态 解 见 图 1， 记 前 
切 流 场 和 单 轴 拉 伸 流 场 下 的 稳 态 数值 解 分 别 为 (9,t) M Yell, t) WAIE 


27 
i, ws(0,t)d0 = 1.0 — 3.9 x 1073, (23) 


27 
/ we(O, t)dd = 1.0 — 7.1 x 10738, (24) 
0 
为 检验 计算 结果 的 正确 性 ， 图 2 给 出 了 对 应 迎风 格式 下 的 差分 解 w(9,t) 和 wi(9,t)， 且 有 


27 
/ (0,t)d0 = 1.00424, (25) 
0 


27 
I w.(8, t)dé = 1.013. (26) 
0 


对 比 图 1 和 图 2 可 以 发 现 ，EFG 方法 的 计算 结果 与 差分 方法 的 计算 结果 吻合 .另外 从 (23)- 
(26) 式 可 以 看 出 ， 使 用 EFG 方法 求解 时 ， 概 率 密度 函数 的 归 一 性 可 以 得 到 更 好 的 满足 ， 从 而 
说 明了 本 文 提 出 的 周期 边界 条 件 处 理 方法 不 仅 可 行 而 且 具 有 较 高 的 计算 精度 . 
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0.1 0.2. 
9 6 
(a) 前 切 流 场 下 的 稳 态 解 (b) 拉 伸 流 场 下 的 稳 态 解 
图 1: 无 单元 Galerkin 方法 求 得 的 Fokker-Planck 方程 的 稳 态 解 
0.4 1 
0.8 
0.3 
« g p- 
0.2 a 
0. 
0.1 0.2 
9 6 
(a) 草 切 流 场 下 的 稳 态 解 (b) 拉 伸 流 场 下 的 稳 态 解 
2: ”差分 法 求 得 的 Fokker-Planck 方程 的 稳 态 解 
5.2 二 维 情 形 
考虑 如 下 Poisson 方程 
Use + Uyy 一 一 8r2 sin(2rz) cos(2ry)， (x,y) € (0,1) x (0,1), (27) 
左右 边界 上 给 定 Dirichlet 边界 条 件 
u(0,y) =u(1,y)=0, ye [0,1), (28) 
上 下 边界 上 给 定 周期 边界 条 件 
u(z,1)=u(z,0), Ze [0,1], (29) 
该 问题 具有 解析 解 
uct 一 sin(2rz) cos(27y), (30) 


计算 时 在 求解 区 域 [0, 1]x[0, 1] 内 均匀 地 布置 41x41 个 节点 ， 背 景 网 格 积分 采用 4x4 点 Gauss 积 
分 ， 计 算 结 果 见 图 3. 
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(a) 数值 解 等 值 线 (b) 解析 解 等 值 线 





(c) 数值 解 (d) 解析 解 


图 3: 二 维 Poisson 方程 求解 结果 


若 记 数 值 解 为 迪 ， 并 基于 L 范 数 定义 如 下 误差 
1 1 t 
e = luexact _ yh le = uexact _ ah) 2 gad 2 i 
l lza = [f f woe - uaran] (31) 


计算 所 得 误差 为 s = 0.00184641. 计算 结果 表明 本 文 提出 的 周期 边界 条 件 处 理 方法 对 二 维 问题 
仍然 有 效 . 


6 ”结束 语 


本 文 针对 具有 周期 边界 条 件 的 微分 方程 ， 提 出 了 EFG 方法 中 周期 边界 条 件 处 理 的 Lagrange 
RFE. 该 方法 将 周期 边界 条 件 视 为 约束 ， 写 出 与 之 等 价 的 积分 形式 ， 然 后 使 用 Lagrange R 
子 法 将 其 引入 到 Galerkin 变 分 方程 中 ， 从 而 使 周期 边界 条 件 在 积分 平均 意义 下 得 到 满足 ， 对 一 
维 Fokker-Planck 方程 和 二 维 Poisson 方程 的 求解 ， 验 证 了 该 方法 的 可 行 性 . 并 且 数 值 计 算 结果 
表明 : 该 方法 具有 较 高 的 计算 精度 ， 在 周期 性 问题 的 求解 中 保留 了 无 单元 Galerkin 方法 高 精度 
的 优点 . 
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Imposing Periodic Boundary Conditions in Element 
Free Galerkin Methods 
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Abstract: In this paper, the Lagrange multiplier method has been used to impose periodic boundary 
conditions in the element free Galerkin method. Numerical results indicate that the method maintains 
the high accuracy property of the element free Galerkin method in periodic problems. In addition, 
based on the similarity between the method used in this paper for periodic boundary conditions and 
the Lagrange multiplier method for essential boundary conditions, it is convenient to impose periodic- 
essential mixed boundary conditions. 
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